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ABSTRACT 


Two independent trajectory software systems are used to perform the TOPEX/Poseidon operational orbit 
etermination and propagation: the Goddard Trajectory Determination System (GTDS) at the GSFC/FDF 

°^ bl t Pr ! CiSi ° n Tra, ' ect0ry System (DPTRAJ ) at J PL - GTDS * used for operational tracking and 
TDRS-based orbit deterimnation. DPTRAJ is used for ephemeris generation necessary to conduct day-to- 
day mission operations This paper describes the DPTRAJ /GTDS trajectory comparison analysis 

conducted jointly by JPL and GSFC to ensure the compatibility of these two independent trajectory 
software systems. } 7 


INTRODUCTION 


The Ocean Topography Experiment (TOPEX) /Poseidon spacecraft was launched on August 10, 1992 to 
study ocean circulation and its interaction with the atmosphere, to improve our knowledge of climate 
change and heat transport in the ocean, and to study the marine gravity field. These objectives are 
accomplished through accurate mapping of the ocean surface with a dual-frequency on-board altimeter 
and precision orbit determination. 

Two independent Orbit Determination (OD) processes are associated with the mission. A Precision Orbit 
Determination (POD) process which is used to support analysis of the altimeter data, and an Operational 
Orbit Determination (OOD) process which is used to support the daily satellite operations. This paper is 
concerned only with the utilization of the OOD solutions in daily operational navigation The OOD is the 
responsibility of the Goddard Space Flight Center (GSFC) Flight Dynamics Facility (FDF). Using tracking 
data from the Tracking and Data Relay Satellite (TDRS) System (TDRSS), the FDF produces 
TOPEX/Poseidon and TDRS state vectors for transmission in the Extended Precision Vector (EPV) 
message format. These EPV solution sets are transferred to the Jet Prbpulsion Laboratory (JPL) via 
National Aeronautics and Space Administration Communications Network (Nascom) to be used by the 
Navigation Team (NAVT) as initial conditions for propagating the Operational Orbit Ephemeris (OOE). 
Operational navigation support procedures have been developed to ensure the compatibility of the FDF- 
estimated TOPEX/Poseidon and TDRS state vectors and the NAVT-generated OOEs (Fig. 1). The 
objective of this paper is to present the results of this activity for only TOPEX/Poseidon trajectories. 

Two independent trajectory software systems are used to perform the above task: the Goddard 

Trajectory Determination System (GTDS) at the GSFC/FDF and the Double Precision Trajectory System 
(DPTRAJ) at JPL. GTDS is used for operational tracking and TDRS-based OD. DPTRAJ is used for OOE 
generation necessary to conduct day-to-day mission operations. 
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Figure 1. OOE Generation & EPV Validation Procedure 


This paper describes the GTDS/DPTRAJ trajectory comparison analysis conducted jointly by the authors 
and dieir colleagues The paper discusses both the TOPEX/Poseidon mission requirements as they 
«1^ to OO^arcuracy requirements and the error budget developed to meet these requirements. The 
operating mocedur^useo^to generate orbit solutions at GSFC/FDF, to transfer ^ solutions to JPL, 
and to process the solutions at JPL are also addressed. A description of each force model enhancement to 
GTDS Ld DPTRAJ motivated by the TOPEX/Poseidon mission is presented followed by a discu^ion of 
the DPTRAJ/GTDS comparison tests and test results. Finally, the lessons leame rom J 
GSFC/FDF experiences; providing joint flight dynamics operational navigation support for 
TOPEX/Poseidon are presented. 


^rriTPACY REOl 
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ON THE OOE 


The TOPEX/Poseidon project has imposed several accuracy requirements on TOPEX/Poseidon 
operational navigation support. The primary driver behind these requirements is a ±1 km error tolerance 
on the equator crossings of the satellite ground track to maintain ground track repeatability. Orbit 
Maintenance Maneuvers (OMM) used to maintain this ground track must be planned and evaluated to a 
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commensurate accuracy level based on the OOE's. More specifically, a 30-day OOE must have a la error 
of no more than 250 meters in equator crossing location. To ensure this level of accuracy, an error budget 
was prescri ed (Ref. 1) that apportioned the overall error allowance for the OOE’s among the identified 
error sources. The one error source of interest here is the trajectory software modeling errors. Ref (1) 
mdKates that the la value of this error must be no more than 40 meters in equatorial crossing longitude 

r7 S ' ■ A 5 me ' K (OUl °' * he 30 120 '“"Situde difference ». equafor croLgt 

after 30 days was allocated as a derived requirement for JPL-FDF trajectory prediction software 
comparisons. This 5 meter tolerance was chosen during software development as an achievable objective. 

is more of a goal than a requirement. Joint mission support by the NAVT and GSFC/FDF thus 
demands a high level of consistency between GTDS and DPTRAJ. Both systems must, for example, 
model atmospheric drag and geopotential forces consistently and utilize common drag and solar 
radiation pressure spacecraft cross-sectional area profiles. 


Table 1. Predicted Ephemeris Error Budget 


ERROR COMPONENT Equator Crossing Error at 30 days 

(1 Sigma Random or max. systematic) 
(m) 


Definitive OD (all sources) 

75 

Prediction error (Nature's 

unpredictability after the 
definitive OD interval: density, UTI) 

130 

Prediction Trajectory Software 
Modeling Errors 

40 

Maneuver Execution 

70 

Geopotential tuning limitations 

10 

Total error (Uncorrelated errors) 

171 

Total error 

221 

Allowable error 

250 (TBR1 


Margin available: (Uncorrellated) 
(Correlation - 1) 


BACKGROI IMP 

^dGTDS ef ^ r !, W ff S Hit" 6 durin 8 t ^f mission development phase to ensure the compatibility of DPTRAJ 
and GTDS. This effort began in 1987 (Ref. 2 and 3) by identifying the force models to be used in the two 

SS IrwToTird E fb /POSei d d °, n ° Pe “ Over the yearn man^ cat havtbet 
siaDiisnea to allow for a model-by-model comparison between DPTRAJ and GTDS. 
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Although extensive effort was made during the software development phase to eliminate ^consistencie 
be ween the two systems (Ref. 4 and 5), it was not feasible to eliminate all inconsistency. DPTRAJ and 
GTDS utilize different modeling design implementations, input/output interfaces, numerical strategies, 
3^. ^m“oaSdifferen. force models. For example, DPTRAJ uses a variable-step integrator wrth *e 
tertsTrial dynamical time as an independent variable whereas GTDS uses a t.xed-st. <J ” f 

time as an independent variable. In addition, DPTRAJ uses a conical umbra shadow model for so a 
radiation pressure modeling whereas GTDS uses a cylindrical umbra shadow model. A paramount 
objective of this work was to quantify the effect on OOE error of system inconsistency and to -ensure 
conformity with the error budget. To this end, ten tests were designe^ J*rfo^du^ thep 
mission phase to allow for a model-by-model comparison between DPTRAJ and GTDS. Alter la 
upon .L! dfiolery of an unanticipated thrust-like perturbation, an eleven* test was devtsed to ensure 

consistency of thrust modeling. 

QOT T TTTONS AT FDE 

TOPEX/Poseidon operational navigation is supported by TDRSS which consuls of a set c >ffive 
geostationary relay spacecraft, called TDRSs, the White Sands Ground Terminal if RouSy 

Sands New Mexico and the Bilateration Ranging Transponder System (BRTS). Routinely 
TOP^/Pr e idon is supported by two of the five TDRSs, an "east" and a "west" TORS, which provide 

and two-way Doppler tracking measurements, as well as spacecraft 
commanding^ telemetry support: BRTS provides ground-based range and two-way Dopp 
measurements for use in TDRS orbit determination. 

The FDF generates OOD solutions for TOPEX/Poseidon using the GTDS. GTDS uses a batch weighted 
Sst squares estimation algorithm, in conjunction with TDRSS one-way and two-way Doppler tracking 

Sglfalo^trlck tSTcah^g JSSSTtoS ^iSutions are geneJated prior to ^^^^rately 
from the TOPEX/Poseidon OOD solutions using BRTS tracking measurements. The TDRS OOD 
solutions are then used in generating the TOPEX/Poseidon OOD solutions. 

DATA TRANSFER TO IEL 

. ff fmm C^VC /FDF to 1PL via Nascom in series of 4800-bit data blocks 

EPV state vectors are transmitted from 10 ji l via T rsnuv /p AC pirlnn OOD 

SSSSSaSILSi: 

TOPFY NAVT at TPL These vectors have epochs at the start and end > 

TOPEX N A V 1 at j i l. nesevt ? , Wednesdays, the vector with epoch 

=E 

as well as to quality assure the OOEs generated by the NAVT. 

designed to raise the semi-major axis to rf ■*•*« vectors 

band. In support of these in addition to .he routine OOD. 

OOD arc. These state vectors are used to assist me inav i in evmum h r 

calibrating the thrusters. 
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IFL PROCFSqpqr; 

An ephemeris file for either TOPFY /Pr^j 

state vectors provided by GSFC/FDF. In' addition to EPV ® ei ? erated by DPTR AJ for each set of EPV 
file for DPTRAJ d ml ° Pr0CCSS by the soft,vare modu le GfNDR^whic^rod^cet^ a NamXw^e 

over the required time span. Per validation, the emaiS," EPV,!® Sa ? Uite ' S ei 5“ ali °" 

state vectors extracted from the ephemeris. 8 EPV compared with their corresponding 


PROPA GATION Monpr q 


^ Pore, don p^els used 
precision prediction can be found in Ref. 6): " maneuver requirements for subsequent 

• Geopotential Model 

STuCluy!! Goddard Model (GEM)- 

degree and order 70. A .3L, 2„ x 20 "? ^"-wilh respeeHp a body-fixed frame up 
computational limitation at JPL. S1 ° n 1S USe ° r °P eratlonal navigation because of 

• Lum-Solar Gravity 

perhfrbing C “ ^ m ° deled ad «P»tely b Y considering these 

• Solid Earth Tides Model 

compensate for the defomatfon of^eTohd porton^ ° f * he ^opotenUal model to 

the Sun and the Moon. The model includes a^ae angle hpJ! by the combuied tidal effects of 

of the disturbing body and the stretching axis Th e model alf^ be az,muthal component of the position 
proportionality constant for the effect § As 3 W nUmber which ser ves as a 

additive adjustment to the gravitational fon^ OT^tlw^pacecraft ^ DPTRAJ ' the m ° del yields - 

# Atmospheric Drag 

major axis reduction. the ° rbital decay m terms of semi- 

relative to the atmosphere. Density is a conflicated function of ° d , enSlty and the velocity of the satellite 
geometric parameters, and diurnal, annual, and latitudinal seaso^l^ 311(1 g eo m a gnetic activity, satellite 
use the same solar and geomagnetic activity data suonbeTh T ^f nataons - Both Dp TRAJ and GTDS 
Admmistration (NOAA). The Jacchia-Roherls ahnosph^ric 
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• Solar Radiation Pressure: 

„ , oT nPFX /Poseidon that exceed those of atmospheric drag, 

The solar radiation pressure (SRP) ^as effec so ° 4eU . The effect of the numerical integration due 
however, this perturbation can be modeled ^ " wrbatio „ when the satellite passes through 

,o the extremely rapid changes m ^ rad,a '‘°" n S model that allows for no integrator restarts has been 
S^n^eSi^soLare. GTDS does no. restart the integration either upon entry 

exit from its cylinder shadow model. 

• Variable Mean Area Model 

The variable mean area (VMA) model allows for a > ™ a ^”3 for , totinc. 

purpose of computing per.urbat.ons due drive „ „ y . par ame.er called (S'. which is 

SRP and atmospheric drag area pro . P ^ spacecraft orbital angular momentum 

the compliment of the angle between the ^-sun ^ ^ ^ sleering yaw -, , table of 

profile is overridden with constant area values. 

• Thrusting Effects 

. ,. . . nrhital decaY levels about 60 times larger than could be 

Shortly after launch OD solutions . md.cated o b ^ a presence of body-fixed 

explained by atmospheric drag (Ref. 7 )- ' caused either orbita l decay or boost dependmg on 

residual along-track forces comparable to d g Consequen tly, plans with the FDF were made to 

m'cTa continuous finite bun, with duration e q ua, to the length o, 
the OD arc and force equal to (1 + t) micro Newton. 


tp AJECIQR Y rnMPARlSQNJlSIS 


| K A 1 ELV- i yixi v- > m 

In all eleven traiectory comparison tests were conducted in prepamtion <* 

Lugh no. identical to the tests originally perfo ^^“"^operational support configuration 
TOPEX/Poseidon. The current set o tes s is in igina i tests were performed). The tests provi e a 

(the configuration has evolved Each teSt * 

model-by model comparison between GSFC ZnZ ^ih DPTRAJ and GTDS. A comparison of the 

- —s «■ “ amined ai resu,ar m,e 

during the 30-day period. 

operational support configuration. 

The tests were as follows: 
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1. Point mass earth 


Ephemeris 


2, Earth gravitational perturbations 


Z * » « .» Additional e P he m e ri s 
fixed coordinate transformations. P implementation and in inertial to earth- 


3. 


Earth, Sim, and Moon gravitational perturbations 


confirm consistency in soll^^l^r eph^el^s. ^ 3X1(1 M °° n ' agreement for ^ test w °tdd 

4. Solid-Earth Tides 


This test adds to Test 3 the solid-Earth tides effects on the geopotential 
iscrepancy would be negligible because of the simplicity of the tides model. 


Additional ephemeris 


5. Expanded Gravity Field 

wmewhat nmre*^^ that^aa io rxii ? UP *° 26 X ^ EphemeriS dfaa W 

6. Solar Radiation Pressure 

This test adds to Test 1 perturbations due to SRP on the satellite AHdiHrmai Q u . ,. 

would arise from shadow model differences (conical umbra for DPTOaI V ' ph , em f ris discrepancy 

smaU difference is expected in this test. computations. As Ref. (3) indicates, a 


7. 


Variable Mean Solar Radiation Pressure Area Model 


This test adds to Test 6 the VMA model for thp cpp * jj*.- i , . 

negligible because VMA implementations in DPTRAJ and in W ° Uld ^ 

8. Atmospheric Drag 

This test adds to Test 1 perturbations due to atmospheric drag The Tacchia R n Wc ^ , . , 

cons.anlspacUc^^^ differences. ibis L, 


9. Variable Mean Atmospheric Drag Area Model 


would be negligtbuTbecause X 
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10. Combined 

prescribed by the error budget. 

11. Operational Test 

“ rol^s?^s t,o rr d n, ! y o rdun, t ? h " 

perturbations was not identified until after launch. 
flf ST TT ,TS AND FUT T TFF. CQNSIPERAJIQNS 

Excellent m odel-b y -mode. agreetnen. between GTDS «d ™J ^^*tl^”a7*-ted 
system to be used to, good orbit prediction 

longitude at equator crossmgs af y above figure is not necessarily a worst case. 

SSSS two systems exceeded the 5 meter goal. Figure (3) shows a 

one-year statistics of the longitude difference. 


Case 11 (Case 10+Thrust) 

I 

Case 10 (Combined Cases 4+7+9) 

Case 9 (Case 1+Drag w/VMA) 

Case 8 (Case 1+Drag w/Sphere) 

Case 7 (Case 1+Solar Pressure w/VMA) 

Case 6 (Case 1+Solar Pressure w/Sphere) 

Case 5 (Case 4+26x26 Gravity) 

H 

Case 4 (Case 3+Solid Tides) 

Case 3 (Case 2+Luni-Solar) 

Case 2 (Case 1+20x20 Gravity) 

Case 1 (2 Body) 
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Figure 2 Trajectory Comparison Results 
JPL(DPTRAJ)-GSFC(GTDS) 
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Figure 3. One-Year Trajectory Comparison Results 


LESSONS T.F A RMrp 

Below are the lessons learned while performing the DPTRAJ/GTDS trajectory comparisons: 

Lesson 1 : Perform the trajectory calibration analysis as early as possible. 

Calibrating any two complex software systems such as DPTRAJ and GTDS can be a tedious and time 

Pr ° CeSS b6CaUSe ° f 1116 lar S e number of variables invoLd. Recognizings fact the 
the NAVT mitia.ed .he trajectory ealibrahoa effort almost thtee yearn befoTia^fch *" Ms 

betw^D^T^^ 0 " lha " SUffi “ n ' ,ime 10 iden " y - a " al> '“ and discrepancies 


Lesson 2: 


When software is developed, all constants should be user modifiable. 


When performing a trajectory calibration analysis, it is critically important to ensure that both software 
ystems use the same modeling constants. It is equally important that the user be able to easily = 
any constants which must be changed for compatibility. In the case of the DPTRAJ/GTDS trajectorJ 

2,d not' h n a' H ? nSt T tS Were easiJ y modified ' s ™ e constants where input by the usZ 

and not hardcoded with.n the software. Had the constants been hardcoded, the trajectory calibration 
analysis would have required significantly more time and effort. Y CaIlbrahon 


177 









acknowledgments 

The research described in this paper was carried out by the Jet Propulsion Laboratory, California Institute 

rf jXo^y unde, a contract with the National Aeronautics and Space Admm.s,rat.on The author 

would like fo acknowledge Ray Frauenholz for his valuable discuss.ons and suggest, ons and Bryan 

Brown for his useful input to and review of this paper. 

REFERENCES 

1. TOPEX/Poseidon Mission and Systems Requirements, 633-103, April 1989 (Internal Document). 

2 Frauenholz, R.B., "Plan and Schedule for Trajectory Propagation Software Comparison Between 
GSFC and JPL for TOPEX", IOM 314.5-1061, 16 January 1987, JPL (Internal Document). 

3. Frauenholz, R.B., "Summary of GSFC/JPL Trajectory Propagation Software Comparison for 
TOPEX", IOM 314.5-1061, 2 July 1987, JPL (Internal Document). 

4 Guinn J., A Presentation at TOPEX/Poseidon Navigation Development Review 4, April 10, 1991 
(Work is done with Mike Nemesure's help), Internal Document. 

c c alama A . A Presentation at TOPEX/Poseidon Operational Navigation Meeting No. 9 at GSFC, 
November 19-20, 1991 (work is done with Mike Nemesure's help). Internal Document. 

6 Bhat, R.S., R.B. Frauenholz, and P.E. Cannell, "TOPEX/Poseidon Orbit Maintenance Maneuver 
Design". AAS/AIAA astrodynamics Specialist Conference, Stowe, Vermont, Aug. 1989. 

7 Fraueholz R B T W, Hamilton, B E. Shapiro, and R.S. Bhat, The Role of Anomalous Satellite-Fixed 
Accetahons L TOPEX/Poseidon Orbit Maintenance", AAS/AIAA astrodynam.es Specahst 
Conference, Victoria, B.C., Canada, Aug. 1993. 


178 



